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Abstract. We consider a spin-1/2 Bose-Einstein condensate with equal 
Rashba and Dresselhaus spin-orbit coupling. After reviewing some rel¬ 
evant features of the quantum phases of the system, we present a short 
study on how their properties are changed by the presence of non-zero 
magnetic detunings and spin-asymmetric interactions. At small values 
of the Raman coupling and of the magnetic field the so-called stripe 
phase occurs, which displays both superfluidity and periodic density 
modulations, in analogy with supersolids. We finally review a recent 
proposal (Phys. Rev. A 90 , 041604) to improve the visibility of the 
fringes, based on the space separation of the two spin components into 
a 2D bi-layer configuration and on the application of a 7r/2 Bragg pulse, 
and we show that this new configuration also yields a sizable increase 
of the stability of the stripe phase against magnetic fluctuations. 


1 Introduction 

The recent experimental realization of artificial gauge fields on neutral atomic Bose- 
Einstein condensates (BEC) [ 1121314 ] represents one of the most important achieve¬ 
ments in the physics of ultracold gases. In the last few years, the nontrivial properties 
of such systems have attracted a broad interest, resulting in a wide number of papers 
devoted to this subject (see, for example, the reviews 060 and references therein). 

A very relevant class of quantum gases coupled to synthetic gauge fields is rep¬ 
resented by spin-orbit-coupled configurations. The first experimental implementation 
of spin-orbit coupling on a neutral atomic gas was performed by the NIST team in 
[D, where they managed to realize a BEC with equal Rashba [5] and Dresselhaus 
[9j spin-orbit couplings. The phase diagram of this system exhibits novel quantum 
phases, which include a stripe phase and a spin-polarized plane-wave phase mm- 
The stripe phase is characterized by periodic modulations of the density profile re¬ 
sulting from the spontaneous breaking of translational symmetry, similar to what 
happens in supersolids [12] . Although experiments have already been made in the 
relevant range of parameters, a direct evidence of such modulations is still lacking, 
mainly due to the smallness of the amplitude and of the wavelength of the fringes. 

In the first part of this paper we review the ground-state properties of a spin-orbit- 
coupled BEC in uniform matter, and we study how such properties are affected by 
the introduction of a non-zero magnetic detuning and of spin-asymmetric interaction 
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strengths. We then discuss a combined procedure to make the stripes visible and 
stable BE thus allowing for a direct experimental detection. 


2 The model and the quantum phases 

2.1 Single-particle Hamiltonian 

We consider a spin-1/2 BEC with the kind of spin-orbit-coupling first realized by the 
NIST team in [3] ■ When written in a locally spin-rotated frame [T4| , the single-particle 
Hamiltonian describing the system is given by (we set Ti = m = 1) 

1 2 

h sp = ^\(Px ~ k 0 cr z ) + 

It accounts for the application of two counterpropagating and polarized laser fields, 
with wave vector difference ko chosen along the x direction, in the presence of a 
nonlinear Zeeman field. The two lasers provide Raman transitions between the two 
spin states, with Raman coupling strength given by 17. The effective magnetic field S 
is given by the sum of the true external magnetic field and of the frequency detuning 
between the two lasers (see, for example, !14j). The spin matrices entering the single¬ 
particle Hamiltonian |l|) are the usual 2x2 Pauli matrices. It is worth pointing out 
that the operator p entering ([I]) is the canonical momentum — iV, with the physical 
velocity being given by v± = p =p koe x for the spin-up and spin-down particles. In 
terms of p the eigenvalues of ([I]) are given by 

* ± <P ) = ^ + (2) 

where E r = kZ/2 is the recoil energy. The double-branch structure exhibited by 
the dispersion pi) reflects the spinor nature of the system. A peculiar feature of the 
dispersion pj, when S = 0, is that it displays, as a function of p x , two degenerate 
minima at ±fco\/l ~ (H/4E r ) 2 , both capable to host Bose-Einstein condensation. 
Notice that the wave vectors corresponding to such minima differ from ±k 3 if 17 ^ 0, 
and vanish for 17 = 4 E r . For larger values of 17 the gas is in the single-minimum 
phase, where all the atoms occupy the p = 0 single-particle state. 


2 , n 

p±\ + j a * 


(i) 


2.2 Many-body ground state 

The peculiar features of the single-particle dispersion p| are at the origin of new 
interesting phases in the many-body ground state of the BEC. For a gas of N particles 
enclosed in a volume V, in the presence of two-body interactions, the mean-field 
interaction Hamiltonian takes the form 


Hint = / d-' 


^yM r ) 2 + ^Y n i ( r ) 2 + m r R( r W( r ) 


( 3 ) 


where g aa ' = ATTa aa ' (a, a' = t, 4-) are the coupling constants in the different spin 
channels, fixed by the corresponding scattering lengths a aa ', while n-j-^ are the densi¬ 
ties of the two spin components. The quantum phases predicted by mean-field theory 
depend on the value of the relevant parameters ko, 17, S and the interaction parame¬ 
ters G 1 =n (g n + g u + 2 g n ) /8, G 2 = n (g n + g u - 2g n ) /8, G 3 = n (g n - g u ) /4 
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Fig. 1 . Phase diagram of a spin-orbit-coupled BEC. The color represents the value of fci/feo- 
The white solid lines identify the phase transitions. The quantity fv c ' = E r /(g'y) is the 
density at the tricritical point where the three phases meet. The diagram corresponds to a 
configuration with 7 = (g — g-fi)/(g + < 7 - 1 - 4 .) = 0.0012 consistent with the value of [ 2 ]. 


eh, with n = N/V the average density. In uniform matter, the ground-state wave 
function can be determined through a variational procedure based on the ansatz 


miiMhJ 


<F(r) = \fn 



( cos 9 + \ 
sin 9 + J 


e ik + x + C. 


f sin 0 _ 
cos 9- 


e 


—ik—x 


(4) 


where C + and C_ are coefficients satisfying the normalization constraint | | 2 + 

|C _| 2 = 1, and k± represent the momenta at which Bose-Einstein condensation takes 
place. For a given value of ko, f2, 6 and the two-body interaction strengths, the values 
of the variational parameters C±, k± and 9± can be calculated through a procedure 
of energy minimization, including both the single-particle 0 and the interaction © 
terms in the Hamiltonian. In particular, energy minimization with respect to k± yields 
the general relationship 2 9± = arccos (k±/ko) fixed by the single-particle Hamiltonian 
0 . Once all the variational parameters have been determined, one can calculate key 
physical quantities like, for example, the momentum distribution, accounted for by 
the parameters k± , the densities n | and n 4, of the two spin components, the total 
density n = = n-f + nj., the spin densities Sk = ^ with k = x, y, z and the 

corresponding spin polarizations (<jfc) = IV -1 f d 3 r s*,. 

The full variational calculation has been performed in EH, where the case of 
spin-symmetric coupling constants = g ^ = g and zero detuning 5 was mainly 
considered. In this situation one finds k .|_ = fc_ = k\ and 9+ = 9_ = 9 for symmetry 
reasons. The ground state was found to be compatible with three distinct quantum 
phases; the corresponding phase diagram is shown in Fig. [T] 

(I) Stripe phase. For small values of the Raman coupling fi and g > g-jq, 
the ground state is a coherent superposition of the two plane-wave states e ±lklX 
with equal weights (|C+| = |C_| = l/\/2), yielding a vanishing longitudinal spin 
polarization (a z ). The most striking feature of this phase is the appearance of density 
modulations in the form of stripes according to the law 

cos {2k\X + (j)) , 


n( r) = n 


2 (2E r + Gi) 


(5) 
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with the periodicity of the fringes n/ki fixed by the wave vector k\ = k^^/l — [17/(2(2 E r + Gi))] 2 . 
These modulations appear as the result of a mechanism of spontaneous breaking of 
translational invariance, with the actual position of the fringes being given by the 
value of the phase (j>. The contrast in n(r) is given by 


^max 

^max T ^min 


n 

2(2E r + G\) 


( 6 ) 


and vanishes as 17 —► 0 as a consequence of the orthogonality of the two spin states 
entering Eq. (4| (in this limit 9 —> 0 and k\ —> ko). It is worth mentioning that the 
ansatz, Eq. (14), for the stripe phase provides only a first approximation, which ignores 
higher-order harmonics caused by the nonlinear interaction terms in the Hamiltonian. 

(II) Plane-wave phase. For larger values of the Raman coupling, the system en¬ 
ters the so-called plane-wave phase (also called the spin-polarized or de-mixed phase), 
where Bose-Einstein condensation takes place in a single plane-wave state with mo¬ 
mentum p = kie x (CL = 0), lying on the x direction (in the following we choose 
k\ > 0). In this phase, the density is uniform and the spin polarization is given 
by (<j z ) = fci/fco wifi 1 &i = koy/1 — [12/(4 (E r — G 2 ))] 2 . An energetically equivalent 
configuration is obtained by considering the BEC in the single-particle state with 
p = —k\e x (G+ = 0). The choice between the two configurations is determined by 
a mechanism of spontaneous symmetry breaking, typical of ferromagnetic configura¬ 
tions. 

(III) Single-minimum phase. At even larger values of 17, the system enters the 
single-minimum phase (also called zero-momentum phase), where the condensate has 
zero momentum ( k\ = 0), the density is uniform, and the average spin polarization 
(cr z ) identically vanishes, while (a x ) = —1. 

The critical values of the Rabi frequencies 17 characterizing the various phase 
transitions can be identified by imposing that the chemical potential g(n) and the 
pressure P = ng(n) — f g{n) d n be equal in the two phases at equilibrium. The 
transition between the stripe and the plane-wave phases has a first-order nature. In 
the low density (or weak coupling) limit, i. e. G\, G 2 E r , the critical value of 
the Raman coupling I7d~ n ) characterizing such transition is given by the density- 
independent expression mm 


< 7 > 

with 7 = G 2 /G i- The transition between the plane-wave and the single-minimum 
phases has instead a second-order nature and takes place at the higher value mi 
I7 (n m ) = 4 (E r — G 2 ), provided that the condition n < is satisfied. For higher 
densities one has instead a first-order transition directly between the stripe and the 
single-minimum phases. We also remark that, if g < g -jq., the stripe phase is ener¬ 
getically unfavorable, and only the plane-wave and the single-minimum phases are 
available. 


2.3 Effects of non-zero detuning and spin-asymmetric interactions 

The results discussed in the previous paragraph can be easily generalized to account 
for the presence of a non-vanishing magnetic detuning S and of spin-asymmetric 
interactions 7 ^ g^. In particular, the effects of an asymmetry in the intraspecies 
coupling constants can be compensated by choosing a magnetic detuning 6 = — 2 G 3 , 
which ensures that the ground-state properties remain the same as in the symmetric 
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Fig. 2. Detuning versus Rabi coupling phase diagram in the experimental conditions of [3]. 
The parameters are E r = 2n x 1.77 kHz, density in the center of the trap no = 1.9 x 10 14 cm ~ 3 
and the scattering lengths given in the main text. 


case m- This holds for any value of G 3 in the plane-wave and the single-minimum 
phases; instead, in the stripe phase exact compensation is possible only if G 3 is small. 

In the most general case of arbitrary 6 and G 3 , the ground-state wave function 
can be still worked out by resorting to the ansatz Q, where now the two wave vectors 
k + and k- can have different values. The resulting phase diagram in the fl-5 plane 
is shown in Fig. [2] It is characterized by the occurrence of a stripe phase, where 
both momentum components of are present (although they can have different 
weights |G+| and |G_|), and by several plane-wave states, having different values of 
the momentum and hence of the magnetization PI- 

The stripe phase occurs only in configurations where > g+\, corresponding 

to the condition of miscibility of the two spin components in the absence of spin-orbit 
and Raman coupling. One can notice that, in systems with almost equal coupling 
constants, it occupies a very small region in the fl-5 plane. For example, in the 
case of the states |f) = \F = 1,to,f = 0) and |4-) = \F = l,m_F = —1) of 87 Rb, where 
apf = 101.41 as and = a-jq = 100.94a#, from Eq. I?]) one finds the value = 

0.19 E r for the critical Rabi coupling, while the critical magnetic detuning needed to 
bring the system from the stripe to the spin-polarized phase turns out to be of the 
order of 10~ 3 £V. The latter is actually proportional to the difference Z\/i between 
the chemical potentials in the two phases. An analytic estimate can be obtained in 
the fl — > 0 limit, where one finds that the stripe phase is favored for values of the 
magnetic detuning 5 such that |i5 + 2 G 3 I < 4 G 2 . At finite values of fl the range of 
values of 5 compatible with the stripe phase is further reduced. As a consequence, a 
tiny magnetic field (arising, for instance, from external fluctuations) can easily bring 
the system into the spin-polarized phases. The stability of the stripe phase can be 
strongly enhanced if one increases significantly the value of G 2 , as we will discuss in 
Sect. 01 

The phase diagram of Fig.[2]also contains three different kinds of plane-wave phase. 
Those on the left region of the diagram, denoted by PW1 and PW2, correspond to the 
regime where the single-particle Hamiltonian has two local minima; PW1 and PW2 
are favored for magnetic detunings <5 larger and smaller than — 2 G 3 , respectively. The 
state PW3 appears instead in the region where the single-particle Hamiltonian has 
one local minimum only. The red dashed lines represent the transition between the 
double-minimum and the single-minimum regimes; in drawing them we have taken 
the presence of spin-dependent interactions into account, yielding some corrections 
with respect to the single-particle results [IB]. 
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3 Experimental perspectives for the stripe phase 


The stripe phase is doubtlessly the most intriguing phase appearing in the phase di¬ 
agram of Sect. [2j It has been the object of several recent theoretical investigations 
rim i ;W)i2 1 . The stripe phase is characterized by the sponta¬ 

neous breaking of two continuous symmetries. The breaking of gauge symmetry yields 
superfluidity, while the breaking of translational invariance is responsible for the oc¬ 
currence of a crystalline structure. The simultaneous presence of these two broken 
symmetries is typical of supersolids [I2l29l30l3i| . It has been shown, among other 
things, to be at the origin of the appearance of two gapless excitations as well as of 
a band structure in the excitation spectrum [33] • 

In the experiments of [3] and m a phase transition has been detected close to the 
theoretical prediction = 0.19 7^ (see Eq. ([t]) ) for the critical Raman coupling 

below which the occurrence of the stripe phase is expected. However, as we already 
mentioned in the introduction, there is still no direct experimental evidence of the 
periodic modulations in the density profile characterizing the stripe phase. The main 
reason is that, in the conditions of current experiments with spin-orbit-coupled 87 Rb 
BECs [4118132) . the contrast and the wavelength of the fringes are too small to be 
revealed. Another issue is represented by the fragility of the stripe phase against 
fluctuations of external magnetic fields, which has already been discussed in Par. |2.3| 
In [13] the authors proposed a procedure to make the experimental detection of the 
fringes a realistic perspective, improving their contrast and their wavelength, and 
increasing the stability of the stripe phase against magnetic fluctuations. 

In order to achieve a larger value of the contrast ([6]), one needs to enlarge the 
range of values of 17 compatible with the existence of the stripe phase. As can be 
seen from Eq. (JtJ) , an efficient way to increase the critical Raman coupling l?^ 11 ) 
is to reduce the value of the interspecies coupling constant g- jq. A possibility is to 
look for hyperfine states characterized by a small (or tunable) interspecies scattering 
length. Here we discuss a different strategy, based on the idea of reducing the effective 
interspecies coupling by means of suitable trapping conditions. In particular, one can 
trap the atomic gas in a 2D configuration, with tight confinement of the spin-up and 
spin-down components around two different positions, displaced by a distance d along 
the z direction. This configuration can be realized with a trapping potential of the 
form 

Kxt ( z ) = y { z ~^) (8) 

produced either through magnetic gradient techniques or via spin-dependent opti¬ 
cal potentials. Assuming a Gaussian profile ip± = (1/ yjnTal)e~^ d ^ 1 2 ^ / 2a * for the 
z dependence of the spin-up and spin-down wave functions, with a z = 1/y/UJz the 
oscillator length along z, the integration over z of the energy functional ([3]) gives rise 
to effective 2D coupling constants g a p given bfl 


fl'tt.H - /7T StbU > 

V 27ra 2 


9n = 


\Z2tto 


-9n e 


-d 2 /2al 


(9) 


In an analogous way one finds that also the effective Raman coupling, to be used 
in 2D, is lowered with respect to the physical coupling 17 according to the law 
17 = e~ d / 4a -17, reflecting the reduction of the overlap between the two wave func¬ 
tions. Hence, the new configuration produced by a tight axial trapping potential with 
a spin-dependent displacement can be described formulating the Hamiltonian in 2D, 

1 In the present section we consider realistic spin-asymmetric interaction strengths ^ 

344 , and we compensate the asymmetry by choosing <5 = — 2 G 3 . 
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Fig. 3. Integrated density profile f dydzn in the striped phase, (a) displays the situation 
without separation of the traps for the two spin components, (b) corresponds instead to 
traps separated along z by a distance d — a z , which increases the contrast of the fringes. 


with the effective Raman coupling given by 17, and the interaction term obtained 
from the functional © with the replacement of the 3D densities with their 2D coun¬ 
terparts f dzn and of the coupling constants with the renormalized values (|9). The 
main advantage with respect to the original 3D problem is that now, due to tire rel¬ 
ative separation of the atomic clouds of the two spin components, the effective g-j-j, is 
reduced with respect to the two intraspecies coupling constants (see Eq. (|9|). As a 
consequence, the ratio 7 appearing in Eq. 0 is now given by the expression 


G2 = fftt + flu - 2 * * m e d2/2a " 
Gi gn + g u + 2g n e~ d2 / 2a i ’ 


( 10 ) 


and is larger with respect to the 3D case, yielding an increased value of the critical 
effective Raman coupling, and thus of the largest reachable contrast of the fringes 
in the stripe phase. For e xam ple, choosing the value d = a z and the 8 'Rb hyperfine 
states mentioned in Par. 2.3 one finds the value 7 = 0.25 for the ratio (10), to be 
compared with the value 7 = 0.0012 for the d = 0 cas<© 

Quantitative predictions for the novel configuration discussed above can be ob¬ 
tained by solving numerically the 3D Gross-Pitaevskii equation. In Fig. [3] we show 
the results for a gas of AT = 4 x 10 4 87 Rb atoms confined by a harmonic potential 


2 Another important change is that, due to the increase of the value of 7 , the critical density 

n^ can be significantly lowered with respect to the value in the d = 0 case, becoming of 

more realistic achievement in future experiments. 









Will be inserted by the editor 



Cl/E r 


Fig. 4. Detuning versus effective Rabi coupling phase diagram in 


the conditions of 


Fig.gn 


with frequencies (cj x ,aj y , oj z ) = 2tt x (25,100,2500) H z, th e scattering lengths a acr > and 
the recoil energy E r equal to those reported in Par. 2.3 and consistent with Ref. [3]. 
Fig. 3 1 corresponds to d = 0, while Fig. [3jo corresponds to d = a z = 0.22 /im. In both 
Fig. 3 1 and Fig. 0 > we have chosen values of the Raman coupling equal to one half the 
critical value to enter the plane-wave phase, in order to ensure more stable conditions 
for the stripe phase. In Fig. this corresponds to 17 = (1/2)17^ n ^( 7 ) = 0.095 E r 
with 7 = 0.0012, while in Fig. ga to 17 = (l/2)e d / 4 “*f 7 ( I-II )( 7 ) = 1.47 E r with 
7 = 0.25. The plotted density corresponds to the ID density as a function of the 
most relevant x variable, obtained by integrating the full 3D density along the y 
and z direction. The figure clearly shows that in the conditions of almost equal cou¬ 
pling constants (Fig. gi) the density modulations are very small, while their effect is 
strongly amplified in Fig. 0 > where the interspecies coupling is reduced with respect 
to the intraspecies values by a factor ~ 0.61. 

The suggested procedure has also the positive effect of making the stripe phase 
more robust against fluctuations of external magnetic fields. Indeed, the reduction of 
the interspecies coupling and the increase of the local 3D density, due to the tight 
axial confinement, yield a significant increase of the energy difference between the 
stripe and the plane-wave phases. For example, in the case considered in the above 
3D Gross-Pitaevskii simulation with d = a z (Fig. §>), a magnetic detuning of the 
order of 0.3 E r is needed to bring the system into the spin-polarized phase (see the 
diagram in Fig. g, while in the abse nce of displacement (Fig. g) the critical value 
is much smaller (~ 0.001 E r , see Par. 2.3 and Fig. g 

Let us finally address the problem of the small spatial separation of the fringes, 
given by ir/ki, which turns out to be of the order of a fraction of a micron in standard 
conditions. One possibility to increase the wavelength of the stripes is to lower the 
value of k 0 by using lasers with a smaller relative incident angle. In the following we 
discuss a more drastic procedure which consists of producing, after the realization 
of the stripe phase, a 7t/2 Bragg pulse with a short time duration (smaller than the 
time l/E r fixed by the recoil energy), followed by the sudden release of the trap. This 
pulse can transfer to the condensate a momentum ks or — ks along the x direction, 
where ks is chosen equal to 2fci — e with e small compared to k\. The 7r/2 pulse has 
the effect of splitting the condensate into various pieces, with different momenta. The 
situation is schematically shown in Fig. [5] for the spin-down component, where the 
initial condensate wave function, which in the stripe phase is a linear combination 
with canonical momenta ±fci, corresponding to momenta ko — k\ and ko + k\ in the 
laboratory frame, after the Bragg pulse will be decomposed into six pieces. Two of 
them, those labeled in the lower part of the figure with momentum ~ 0, will be 
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0 ~ +2fco 



Fig. 5. Schematical description of the splitting of the spin-down component of the stripe 
wave function into different momentum components caused by a 7r/2 Bragg pulse transferring 
momentum 2k\ — e. 



-20 -10 0 10 20 


x [/im] 



Fig. 6. Integrated density profiles f dzn (top) and f dydzn (bottom) in the stripe phase, 
in the same conditions as Fig. [3ja, after the application of a n/2 Bragg pulse transferring 
momentum ±1.8hki. 


practically at rest after the pulse and are able to interfere with fringes of wavelength 
27r/e, which can easily become large and visible in situ. It is worth noticing that 
these two latter pieces originate from the two different momentum components of the 
order parameter Q in the stripe phase and involve 1/3 of the total number of atoms. 
The corresponding interference effect would be consequently absent in the plane-wave 
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phase, where only one momentum component characterizes the order parameter. The 
other pieces produced by the Bragg pulse carry much higher momenta and will fly 
away rapidly after the release of the trap and of the laser fields. In Fig. [6] we show 
a typical behavior of the density profile obtained by modifying the condensate wave 
function in momentum space according to the prescription discussed above. 


4 Conclusions 

In this paper we have reviewed the ground-state properties of a spin-orbit-coupled 
BEC, focusing on the effects of the presence of a finite magnetic detuning and of spin- 
asymmetric coupling constants. The phase diagram includes a stripe phase, which is 
characterized by the presence of periodic modulations in the density profile. In the last 
part of the paper we have discussed a combined procedure to increase the visibility of 
such modulations and the stability of the stripe phase against magnetic fluctuations, 
thus favoring the exploration of this intriguing configuration in realistic experimental 
conditions. 


The author acknowledges stimulating discussions with Yun Li and Sandro Stringari. This 
work has been supported by ERC through the QGBE grant and by Provincia Autonoma di 
Trento. 
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